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l.INTRODUCTION 

The vorticity equation for an ideal fluid on a two-dimensional manifold Ai is 

C + {C,^} = o, (1) 

where ( is the vorticity and t(j the stream function related to ( by 

C = -A2V'- (2) 

A2 is the Laplace operator on A^ and {/, (7} is the Poisson bracket of / and g. 

This equation has, of course, been the subject of numerous studies over the 
years. It will be enough to mention the analysis of atmospheric motion (in the zero 
height approximation) and the theory of turbulence. 

A standard approach is to expand t/j in modes of A2 so that (1) becomes a 
coupled mode equation, the coupling coefficients being the structure constants of 
the Poisson algebra, with factors involving the eigenvalues. 

Precisely, define modes, Yq, and eigenvalues, Aq, by 

and expand ( and i/j, 

C = ^ CAaYa = "^ Xci^XaYc,, ^ = J^ tpXaY^, (3) 

a a a 

SO that (1) reads 

XotipXa + XpCXai3^'i/jX(3iljXK = 0, (4) 

where the structure coefficients are defined by 

{Yp,Y^} = Y^CXap^ 

and are given, by orthogonality, as an integral over three harmonics. On the two- 
sphere, ElsasserAl appears to have been the first to write down this integral, al- 
though he does not refer to the Poisson algebra. We will not enter into a detailed 
history of these sphere coefficients. They occur in the work of SilbermanA2 and of 
Baer and PlatzmanA3 in early studies of atmospheric vorticity. It was notedA4, 5, 
sometime later, that the coefficients were proportional to Clebsch-Gordan coeffi- 
cients although the calculation of the reduced matrix element was cumbersome. 

The two-torus, TA2, presents, in some aspects, a simpler situation and its 
Poisson algebra was first discussed by LorenzAG who was concerned to truncate an 
infinite coupled mode system (in the atmosphere) to the simplest nontrivial, finite 
one. The same algebra was later investigated by ArnoldA7, also in connection with 
hydrodynamics. 



The modes on the torus are plane waves, exp(zn.r), (— tt < x < tt, —n < y < n) 
and the eigenvalues are An = nA2, n G ZA2. The expansions of the stream function 
and vorticity are 

V'(r) = y^ipXneXin.r, ({r) = N^CAneAm.r (5) 



with (Xn = An^An. The structure constants are 

C'Arin/ n" = n' A n"5Ann'+n"- (6) 

The problem of mode truncation is a vital one in numerical weather prediction 
and there seem to be no theoretical criteria for its optimum solution. One point is 
that any truncation does violence to the infinite set of conserved quantities (which 
may be taken to be the integrated powers of the vorticity) for equation (1). In 
Lorenz's truncation, for example, only the energy and the enstrophy were conserved, 
and this was considered to be remarkable. 

It is therefore of some interest to develop finite-mode approximations that 
preserve more conserved quantities. Such models are suggested by the factAS, 9 that 
the Poisson algebra structure constants (or, equivalently, the structure constants of 
the area-preserving diffeomorphism group, SDiff(7Vl)), are the limits of the structure 
constants of SU(A^) as A^ tends to infinity, after a simple change of normalisation. 
We can say that the commutator of two elements of the Lie algebra of SU(A^) 
"corresponds", in the limit, to the Poisson bracket of two functions on A1. It can 
be seen that there are at least A^ constants of the motion, corresponding to the 
energy and the A^ — 1 Casimir operators. 

ZeitlinAlO has also suggested and investigated these models in works which 
appeared after our analysis was undertaken. There are certain differences of detail 
and emphasis. A reference to finite models is also made by ZakharovAll. 

2. MATRIX ANALOGUE OF THE VORTICITY EQUATION ON TA2 

Since we are interested in SU(A^) we first present some standard algebraic mate- 
rial regarding its generatorsA12, 8, 9 which we write in the WeylAlS form (see also 
SchwingerA14) 

Jn = ujXnin2/2gXnihXn2 (7) 

where the unitary N x N matrices h and g satisfy hg = ugh and gXN = 1 = 
hXN. We choose A^ odd and u = exp{2ikn/N) where k and A^ are coprime. The 
periodicities, 

^n+iVp = (-l)Ap An+pip2^n, k odd 



and 



Jn+Np = Jn, k even, (8) 



can be used to bring any n onto the N x N lattice, Cn (the unit ceU), defined by 
— {N — l)/2 < Ui < {N — l)/2. We occasionaUy use Coo to denote the entire square 
lattice, ZA2. 

The most popular choice appears to be a; = exp(4z7r/A^) because of the simple 
periodicity (8). 

The Jn satisfy the relations 

JnAf = J-n = (Jn)A-l 
Jn' Jn" = CuAn" A n'/2Jn'+n" (9) 

^Tr( Jn'A t Jn"A) = 5AnV, (n', n" G Cn). 
Splitting up (9) gives the commutation and anticommutation rules, 

2 / kn \ 

[Jn', Jn"] = «T sin f —n" A n' j Jn'+n" (10) 

and 

2 / k7T \ 

[Jn', Jn"]+ = T COS f — n" A n' j Jn'+n"- (H) 

Another way of writing (9) is 

JniJn2 = ^eA-Z— — Ai23Jn3 (^ni+nsAns, (12) 

where A123 is the area of the triangle formed by the vectors ni, n2, — (ni + n2). 

If n includes the origin, with Jq = 1, the Jn, n G Cn, form a complete operator 
basis, 

^^JnAtTJ„A = lTrT. 

n 

(SchwingerA14. Incidentally Schwinger chooses k = 1.) 

We noteAS, 9 that, as A^ — *> 00, the structure constants in (10), for finite n' and 
n", tend to those of the torus Poisson algebra, up to a normalisation constant. 

Let V and w be two, time-dependent elements of su(A^). They can be expanded 
in the Lie algebra generators Jn, 

v{t) = ^vXn{t)Jn, w{t) = ^wXn{t)Jn. (13) 

n n 

The summation over n is restricted to the lattice Cn — {O}. If we wish to extend 
the summation to Cat, as we do, vXO is set equal to zero for traceless v. Hermiticity 
is equivalent to the conditions on the coefficients 

f AnA* — vX—n, wXnX* — wX—n. (14) 
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Consider the equation 

v + iN(3[v,w] = (15) 

which we wish to compare to the hydrodynamic equation (1) with v the vorticity 
and w the stream matrices. /3 is a constant that wiU be specified later. There is no 
real signficance to its value since the overall normalisation is actually arbitrary and 
could be absorbed into a redefined time. 

In order to correspond with (2), we require that, as A^ ^ oo, vXn -^ XnwXn. 
Assuming that this has been achieved, the statement about the structure constants 
is that, as A^ tends to infinity, the equation that the coefficients wXn satisfy tends 
to the same equation that the coefficients of the expansion of t(j in torus modes 
satisfy. Then, in the limit, we might hope to identify wXn with t/jXn. This will be 
made more precise later. 

Our intention is to look upon these finite- dimensional models as playing the 
role of consistent Lorenz-type truncations although it is not clear, a priori, whether 
they will prove to be of practical interest. 

We turn first to the relation between v and w, and it is here that we differ from 
ZeitlinAlO. He simply sets vXn equal to XnwXn. We feel that the relation should 
be expressible directly in terms of the Lie algebra elements themselves and it is not 
clear whether this is true for Zeitlin's relation. 

Looking at (2), we require the Lie algebra analogue of the Laplacian. To find 
this we recall the significance of the operators g and /i in (7) as stepping operators 
in the quantum mechanics on the discretised circleAlS, 14, 15. 

It is easy to verify that 

CJn = (^) A2([[/i, Jn],/lA-l] + [[(7, Jn],^A-l]) = -AnJn 

where 

A„(fc) = (^) A2 (.inA2(^) +.i„A2(^)) (16) 

which we recognise as proportional to the eigenvalues of the difference Laplacian on 
the discretised torus {2kn/N)ZN^{2kn/N)ZN. The normalisation factor is chosen 
to give the correct continuum limit. If /c = 1, An ^ An = nA2 as A^ — * cxd for 
fixed n. As a set, the An are independent of k. In fact An(l) = Anp(A;) where the 
components of Up are cyclic permutations of those of n according to 

Hp = kn mod(A^, A^), (17) 

where k is the mod inverse of k, i.e. kk = 1 mod A^. (See e.g. CizekAlG.) 
Another way of expressing this is to say that the discrete ^-function, 

Y^' exp(2zA;7rn.m/A^) 

„^^ (An(^))A. ' 



is independent of k. The prime means that the n = term is to be omitted. We 
note that Jq, the unit matrix, is the zero mode, CJq = 0. 

It might be helpful to remark that the continuum VA2 can be written in terms 
of repeated Poisson brackets, 

VX2'i/j = {{eXix, ifj}, eX—ix} + {{eXiy, ifj}, eX—iy} . 

We take the operator C to be the discrete Lie algebra analogue of A2 = VA2 
so that the generators Jn are the analogues of the modes Yq,, (plane waves on TA2) 
as befits a complete set. The relation between v and w is thus written neatly as 
V = Cw and (15) becomes 

Cw + iNP[Cw,w]=0. (18) 

At this point it is convenient to discuss the conservation properties of (5). We 
first need the fact that C is hermitian, i.e. 

TT{aX]Cb) =Tr {C{aX]b) 

where a and b are elements of su(A^). (Of course aX] = a.) The trace is the finite 
analogue of integration over Ai. 

It is then easy to show that the quantity 

E = —Tr{vw) = 2 5Z KwXnwX-n, (19) 

n 

which we refer to as the energy, is time-independent. 

Also, quite trivially and independently of the relation between v and w, the 
traces of powers of v are conserved 

Si^-T:r{vXl)= J2 vXn^...vXnlCos(J;^A^,„lj, (20) 

ni...n; ^ ^ 

where Ai.j is the area of the (/ + l)-gon in Coo with edges ni, . . . , n^, n^-i-i subject 
to the restriction n^+i = ~-J2i=i '^^^i = (OmodA^, Omod A^). 

If k is even the periodicity of the cosine allows one to replace Ai.j by the area 
of the l-gon, ni, . . . , — J2i ^^ ~ ^^i- {^^ f^ict the whole polygon can be pulled back 
to fit into the unit cell.) 

There are A — 1 independent Si, I = 1,...,A — 1, corresponding to the an- 
ticommutator (11) i.e. to the Casimir invariants constructed from the symmetric 
dXijk SU(A) invariant tensorsA17, 18, 19. 5*1 always vanishes. ^'2 is the enstrophy, 

n. 
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These invariants also arise in the analyses of the generalised Euler equations of 
rigid body motionA9, 20, 21, 22 except that the group there is taken to be SO(A^) 
so that only the even powers remain. 

As A^ tends to infinity, (20) should become the continuum expression, assuming 
that the vXn tend to the (Xn of (5). 

The dynamical equations for the stream element coefficients are 

vXn{t) = AnwXn{t) = - ^ GXun' ,n" K'wXn {t)wXn" (t) (21) 

n',n" 

where the summations are restricted to lie on the lattice Cn — {O} and where the 
coupling coefficients are given by 

GXnn',n" = 2N(3 sin ( -^n A n" j 6Xnn'+n" ■ (22) 

The periodicity (8) has been incorporated by defining the (quasi-)periodic delta, S 
with 

^Ani = y^ 5An + Apj, /seven 

peCoo 

and 

^Ani = ^ {-l)XnAp+pip2dXn + Npi, kodd (23) 

peCoo 

so that n can be restricted to the unit cell. In (23), 1 = n' + n" and the sums are 
actually restricted to p G C3 because adding two elements of Cn can take us only 
to the "nearest neighbour" unit cells. 

As A— i>oo, the GAnn',n" tend to the Poisson algebra structure constants (for 
fixed n' and n") and we expect (21) to turn into Lorenz's torus equation. However 
it is necessary to be careful when taking the A ^ 00 limit. One cannot simply 
substitute the limiting form of (22) directly into (21) because of the behaviour of 
terms for which kn' is of order A, for example. To elucidate this limit we shall 
rewrite the equation in coordinate representation but first, another motivation for 
this particular step will be given. 

The aim is to solve equation (21) numerically for given initial conditions, and 
then to compare with the corresponding discussion of Lorenz, i.e. with a simple 
truncation. Hence there arises a question concerning the appropriate quantity to 
construct once the coefficients have been computed. It is possible to compare the 
coefficients directly but this is sensible only for A large, so making the comparison 
impractical. A stream "function" is needed for a global picture. That is, a contin- 
uous quantity constructed from w{t), for any finite A, that can be compared with 
the conventional stream function ?/;(r, t) after a numerical integration. 

It has been notedA23, 24 that the limit A — > 00 is akin to the transition from 
quantum to classical mechanics with, in these references. An /N playing the role 



of Planck's constant. This suggests that we regard (18) as a Heisenberg equation 
of motion and derive the corresponding classical equation in the standard fashion 
using, say, coherent statesA25. This would give a concrete connection between w 
and the "classical" stream function and will be pursued elsewhere. 

We have not seen a discussion of coherent states in Weyl's finite formulation 
of quantum mechanics although there are several applications of the Wigner phase- 
space techniqueA26, 27, 28 to which we now turn. 

This more formal point of view is provided by the representation of quantum 
mechanics (called "treacherous" by GroenewoldA29) introduced by GroenewoldA30, 
MoyalA31 and others, based on the WignerA32 phase-space distribution, and much 
studied since. The quantum equations are replaced, exactly, by a classical looking 
equation but with the Moyal bracket (actually due to Groenewold) instead of the 
Poisson bracket. 

In this approach, which, a priori, is distinct from the coherent state method, 
one constructs the Weyl- Wigner distributions, Tr(aJn), which are then interpreted 
as the Fourier components of the classical quantity corresponding to the operator, 
a. Usually one starts from a classical quantity and then asks for the corresponding 
quantum operator. This is the well known ordering problem. 

A more general orderingASO, 33, 34 is provided by setting 

a(m) = ^ 5^ Tr (a JAfn) 0(n)eAz-^A;n.m (24) 

n 

(usually k = 1). The Weyl ordering corresponds to O = 1 and then a(m) is called 
the Wigner function. 

If a(r) is to be real when a is hermitian, the function O must satisfy OA* (n) = 
0(— n) and we also want O ^ 1 as A^ ^ oo. Typically, 0(n) is a trigonometric 
function of the product nin2. A gaussian form for O is associated with normal or 
antinormal ordering. 

The quantity that corresponds to the commutator [a, b] is the Moyal bracket 
(if n = 1) and, in this case, the coordinate-space representation of the vorticity 
equation (15) is (the proof is given shortly) 

v{m,t) + {v,w}j^{m,t) =0, (25) 

where the finite Moyal bracket, {a, 6}^ is defined by 

{a,6L(m) = ^ E a(m>(m")sin (^a) , (26) 



m' m 



with 



A= - (m A m' + m' A m" + m" A m) 



the area of the dual triangle with vertices at m, m' and m". 

It is clearly possibleA35 to generalise the Moyal bracket to allow for the more 
general ordering (24) involving O but we will not pursue this point here except to 
say that (7) and (24) show that the different choices for k are related to the ordering 
question. 

As mentioned before, one reason for introducing the coordinate-space repre- 
sentation is that the infinite N limit appears to be more transparent than in the 
mode representation (21) which always remains discrete. We will deduce a value 
for the constant (3. 

Equations (25) and (26) are now derived. The finite Fourier relation we require 
reads 

eXi—rkn.Yn. = >, SXniNp (27) 



A^A2 ^ N ' V 

neCjv peCoo 

and the transform is defined by 

a{m) = y^ aXneXi — kn.ih = N^ aAnpcAz — n.m, 

neCjv neCjv 

1 ■^-^ / \ > 27r , ~ 
aXn = —-—- y a(m)eA— z— -r/cn.m. (28) 

A^A2 ^^ ^ ^ N ^ ^ 

CX*N is the dual lattice. We often identify Cn and CX*n- 

The expressions provide periodic extensions of a(m) and aXn off the corre- 
sponding unit cells. 

As a technical point of some interest, the appearance of the factor of k in (28) 
is related to the use of the eigenvalues An{k) of (16). If we had simply chosen to 
set fc = 1 in (16) (but not in the definition of cu), then the k in (28) must be unity 
too. There is nothing wrong in this, but it would not then be possible to write 
the dynamical equations in purely Lie algebra terms, as we have done in (18). We 
believe this is of more than aesthetic importance. 

As A^ — ^ oo, with r = (27r/A)m and a(m) -^ a(r), the formulae (28) turn into 
a standard Fourier series, 

a(r) = y ^ aAneAin.r, 



\n = - — — — / dra(r)eX—in.r. 



aAn = - — — — / (ira(r) eA— in.r. (29) 

(27r)A2 Jt\2 

and aAn = aAn where n stands for a pair of reordered sets of all the integers. We 
have used the continuum replacement 

^ ' V27r J Jt\2 



It is formally attractive to define the transformed generators, J(ni), by 

J(in) = 2_. JneX—i — kn.m, 



-^^ = Wyy Yl J{m)eXi — kn.m, (30) 



NX2 V ^ 



SO that 

i; ^ > ii\n.T„ = 

NX2 



v = ^vXnJn = jTr^Yv{m)J{-m) 



n 



The J(m) are hermitian, and, as we have defined them, satisfy the relations dual 
to (9) 

TrJ(m) = A^A2, for all m, 

J(m') J(m") = J2 /(m'' in"; m) J(m), (31) 

m 

^T¥(j(m)J(m'))=iVJ^-^,, 

J-mJ{m')Jm = J{m - m). 

These relations hold for all k but an even value would be preferred because of 
the implied simple periodicity of the Jn- 

If k is even, a short calculation using (27) shows that the composition constants 
are given by 

/(m', m"; m) = exp fi^A{ni', m", m)] . (32) 

where A{m' , m", m), = — A(m', m", m), is the area of the triangle (m', m", m) on 
the dual lattice, given before. 

Consider the product of two operators {i.e. Lie algebra elements) 

m' m" m 

all sums being over C^. This defines the *- or Moyal-product. Therefore from (31) 
{a*b){m) = j^ E a(ni')6(m")/(m',m";m). (33) 

m',m" 



Taking the commutator gives 

a. h\ = 

NX2 



1 - 



m 
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where 



[a, 6](m) = {a*b — b * a)(m) = N, a(m')6(m") sin I — — -A(m', m", m) J 

m',m" ^ ^ 



:{a,6L(m). (34) 



A;A3A^ 

Equation (15) can then be written as (25), with (26), as promised, ii (3 — /cA3/27r. 
The continuum hmit of (25) can be checked by replacing m, m', m" by r, r', r" 
respectively, where r = (27r/A^)m etc. (for all k). Then, in the infinite A^ limit, 
m ^ oo, the sums turn into integrals and, just as (28) becomes (29), (26) goes over 
into 

{a,6}^(m) -^ 

^.^ ^A3iVA3 r ^r'rfr"a(r')6(r") sin f^A(r, r', r")l = la, 6}(r), (35) 
TV^oo 47rA5 J \ TV J 

TA2XTA2 

the Poisson bracket, as required, k is assumed to remain fixed and r G TA2. (The 
finite integration ranges could be replaced by infinite ones, in the limit, to give 
precisely the same integrals as in BakerA36.) A(r, r', r") = |(rAr'+r'Ar"+r'''Ar) is 
the area of the coordinate-space triangle (r, r', r"). The conclusion is that equation 
(25) becomes the continuous Euler equation (1). 

The coordinate-space representation also allows us to confirm the form of the 
discrete Laplacian, C The standard expression for the finite- difference Laplacian 
on TA2, scaled to give the correct continuum limit, is 

AA2J(m) = f^) A2 J^ [Ji<^>) " A^)] , 

^ ^ <m> 

where < m > are the nearest neighbours to m. Thus 

AA2J(m) = (^) A2(j(mi,m2-1) + J(mi,m2 + 1) -2J(mi,m2) + 

J(mi + l,m2) + J(mi - l,m2) - 2J(mi,m2)). 
It is easily shown that 

AA2J(m) = f^] A2 ([[/i, J(m)], /lA-l] + [[g, J{in)igX-l]) . 

Thus one can write 

1 1 

Cv = — — V i)(m)AA2J(m) = — — V AA2i)(m) J(m) 
jVA2 ^^ ^ ^ ^ ^ A^A2 ^^ V ; V ; 

m£CjvA* mgCjvA* 
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with AA2i;(m) -^ VA2z;(r), as required. 

The invariants 5"; too can be recast in terms of t'(m). We write 

1 

'^^ " 7iV)A27 5Z v{mi)v{m2)...v{mi)G{mi,...,mi). (36) 

mi...m; 

G is related to a finite Fourier transform, (a Gaussian sum) , by 

1 

G(mi, . . . , m;) = — SymTr (J(mi) . . . J{mi)) , (37) 

where from (30) and (12) 

— Tr(J(mi)...J(mz))= J] eA-i^ I ^ A/n,.m, - Ai., 

ni...n; \ 1 

The symmetrisation on the mi can be performed in various ways. Simply 
reversing their order gives 

^ , , v;^ ^ 2A;7r ^^^ _ /2/c7r . \ , , 

G(mi,...,mi) = 2_^ eA-z— — 2_^A/ni.miCos I — — Ai..j 1 (38) 

ni...n; 1 ^ ^ 

with the mod A^ condition on ^Ui. 

Another formula results on combining the J(m) in (37) using the composition 
law (31) to give 

^TY(J(mi)...J(mO) = 
Y^ /(mi,m2,m3.)/(m3.,m3,m4).../(m(;_i)/,mi_i,mi). (39) 

It is interesting to check that, in the infinite N limit, 5"; becomes the integrated 
power of the vorticity function, that is, up to a normalisation factor. 



Si — > / drv{r)Xl. (40) 

Jt\2 

and it is instructive to carry the limit through completely in coordinate- space after 
the summations over the n^ have been done. 

The behaviour of the function G(mi, . . . , m;) as A^ tends to infinity is required. 
As usual we set r^ = (27r/A^)mi and write G(ri, . . . , r;). Taking the expression (36) 
for the invariant Si, rescaling the mi to the r^ and changing the summations into 
integrations produces 

^i — ' to ^^r,, / dri...driv{ri)...v{ri)G{ri,...,ri). (41) 

[Zn)AZL Jt\2x...T\2 
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The simplest nontrivial example is I equal to three, when the polygons are 
triangles. Then, immediately, from (31) and (32) 

/ 8kn 
G(mi, m2, mg) = NX2 cos ( — — A(mi, m2, ma; 

On rescaling, A(nii, m2, 1113) becomes (A^/27r)A2A(ri, r2, r^) and we can now apply 
another formula given in BakerA36, 

fkN\ , /■ , , , f^kN , ^, , ,, 

hm ] X2u Aoodri . . . dr2,^+i cos (A(ri, r2, rs) + A(ri, r4, rs) + . . . 

iV^oo \ n J J_^ V TT 

. . . + A(ri, r2^, r2^+i)) j F(ri, r2 . . . r2,.+i) 



AoorfriF(ri,...,ri), (42) 

for z/ = 1 to give the desired continuum integral of i'(r)A3 in the infinite A^ limit. 
We note that in (42), F must be a reasonable function. 

The case of any / will now be discussed. Although (39) is completely in coor- 
dinate space it is not in a convenient form for the application of (42). In fact, the 
general form of G(mi, . . . , m;) can be given with no summations. The expression 
depends on whether / is even or odd and, in fact, we shall restrict the discussion to 
odd / for brevity. 

We start from the form (38) and begin by replacing n; by n^ = — J2i -^^ ~ ^^i- 
This is allowed because of the periodicity of the exponential and the cosine. We 
might then just as well rename n^ as n^ and restrict the sum (38) to closed /-gons, 
as mentioned earlier. 

The evaluation proceeds by alternate application of (27) and imposition of the 
resulting d. The choice of which n^ to sum over is crucial to obtaining a simple, 
symmetrical result. It is convenient to first perform a cyclic permutation of the rij 
(under which Aii is invariant) so that n^ becomes n^ where b is the next integer 
after 1/4. Then, performing the sums in the order n^ downwards, we find 

G{mi,. .., mi) = NXl -Icos i^E{{ni,})j , (43) 

where i?({mi}) is given by 
E{{mi}) = ^ A(mi,mi,mj+i) + ^ XI - l{ini - m^+i) A {uij - nij+i) 



i = 2,4,... 
i = 2,4,... 



= A({mJ)+S({mJ). (44) 
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This is a closed form for the Gauss sum. The first summation is the area of (/ — 1)/2 
triangles connected at the vertex mi in a windmill sail pattern. The second sum is 
that of the cross-products of all pairs of vane ends. 

Equation (43) with (44) yields (41) with G given by 

^(ri, . . . , rO = iVA/ - 1 cos (—-{A{{r,}) + S({r,})) j 

.rx, . f t'ikNA. ,2kNB. . ,2kNA. . ,2kNB.\ 
= NXl - 1 cos ( ) cos ( ) - sm ( ) sm ( ) (45) 

\^7r' ^ IT ' ^ TV ^ ^ -K ' ) 

in terms of rescaled quantities. 

A completely immediate application of (42) to (41) is not possible because the 
function F now contains A^. However we note that the effect of the cos {2kN Aj-n^ 
terms in (42) in the N ^ oo limit is to force r^ to equal ri_|_i, (z = 2,4, ...,2/) and, 
also, both to equal ri. Since B has a product structure (see (44)), the cos{2kNB/7r) 
factor can be replaced by unity, in the limit. We also note that the same equation as 
(42), but with a sine (as used in (35) with an extra factor of A^) gives a result of the 
order of 1/A^ and so the second term in (45) goes away. Equation (42) can now be 
applied directly with F(ri, . . . , r;) = i'(ri) . . . v{ri) to give the required continuum 
expression (40). Our discussion of the coordinate representation of the invariants 
ends at this point. 

We can now give a reasonable answer to a question posed earlier regarding the 
appropriate quantity to construct from the matrix w{t) that can be compared with 
the continuum steam function, -(/'(r, t), resulting from a standard mode truncation of 
Euler's equations (1). The Fourier transform (28) suggests the simple interpolation 

w{r,t) = y^ wXnp{t)eXin.r, (46) 

as a possible analogue oiifj^r^t). For convenience, it is this quantity that is plotted 
but it is clear of course that there cannot be a unique quantity corresponding to ifj. 

We note that the Fourier coefficients in (46) are evaluated at the permuted 
points np. This means that, for any k, the plane-wave modes with the smaller 
|ni|, |n2| are associated with the smaller eigenvalues, Anp(A;) = An(l), as occurs in 
the continuous case. A naive application of the N -^ oo limit to the momentum- 
space equations, (21) and (22), does not give the correct result. 

There is a peculiarity in that the coordinate-space treatment of the N ^ oo 
limit is not easily available to us for what appears to be the simplest value of k 
(from the Fourier transform point of view) namely unity. For finite A^, the models 
for different values of the parameter k seem to be distinct. Our treatment shows, 
however, that they will all yield the same continuum limit but we can vouch for our 
coefficients only when k is even. 
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3. MAXIMUM SIMPLIFICATION. 

We now turn to the practical solution of equations (21), along the lines of Lorenz's 
calculationAG. 

Since tj) is real its Fourier coefficients satisfy the condition -i/'AnA* = i/'A— n. 
The obvious finite equivalent is the hermiticity of w, (14). 

LorenzAG notices that if the coefficients are chosen to be real at some initial 
time, they will remain real throughout the time development. The reality condition 
means that wXn. = wX—n and using the symmetries A_n = An and GAnn' n" = 
GA— n_n',-n" = —GXnn",n' we can deduce from (21) that 

d 1 X 

An — {wXn—wX—n) = - \^ {An'—An")GXnn'^n"{'wXn'+wX—n){wXn"—wX—n") 
at Z 

n',n" 

(47) 
showing that if the condition wXn = wX—n is valid for all n at some time, it remains 
valid. 

For A^ = 3 the number of independent real coefficients is four. This is the small- 
est number that we can take for a consistent group theoretical structure. Lorenz 
makes a further identification, that is also propagated in time, and this reduces his 
number to three. In general, the number of coefficients in our maximum simplifica- 
tion is (A^A2 - l)/2. 

For real wXn, w can be rearranged into 



- ^ «;An( Jn + J-n] 



^ 2 

n 

The combination of generators that occurs on the right hand side gives just the 
generators of the subgroup U(SU((A^ + l)/2) (g) SU ((A^ — l)/2) which means that 
the torus has actually been turned into a tetrahedronA37 which might not be un- 
reasonable for discussing atmospheric motion on the whole earth. 

In all our calculations k was set equal to two. The numerical procedure con- 
sisted of choosing an initial distribution of the coefficients wXn and then integrating 
(21) by standard routines. The results for the stream "function" were displayed in 
coordinate space using the Fourier interpolation (46). For each odd value of A^ from 
3 to 31 the results were compared with those for a simple trunctation method using 
the same number of modes. 

The conservation of E and of the 5"; was used as a check of the algorithms and 
algebra. 

The initial configurations are discussed in the next section. 
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4.VORTICES ON A LATTICE. 

For the sake of having something definite, it is interesting to propagate a system 
of lattice vortices. A suitable set of initial stream coefficients for a single vortex 
situated at m = m, would be 



Wi{m) 



neCiv 



exp (2z/c7rn.(m — Yni)/N^ 
AM 



(48) 



They are independent of k. Further, X]m^i('^) ~ ^- "^^^ mode coefficients are 
wXrii = exp (2ik7rn.mi/N)/An{k) if n 7^ and wXOi = 0. For the vorticity, 
vXui = exp (2ik7in.ini /N) if n 7^ with vXOi = 0. In coordinate space, Vi{m) = 
A^A25Amini — 1 showing that the vorticity is mostly concentrated at the point iiij, 
justifying the term "vortex". The smaller, negative value of —1 makes the total 
"integrated" vorticity, X]m'^j('^)' ^^ro, as necessitated by the compactness of the 
domain. However it is not possible to conserve the vorticity located (in some sense) 
at mi, which can leak away. 

As A^ and m tend to infinity, the discrete stream function, Wi{m) of (48), turns 
into the Epstein (^-function, except at a; = 0, y = 0, 



Wi{m) 






exp (in.r) 



nA2 





x/27r y/2Tv 



(2), 



(49) 



in Epstein's notationA38. For simplicity, we have set k = 1 and 111^ = 0. x and 
y are the coordinates of r with x — — 2711712 /N and y = 2TTmi/N. In the limit we 
would regard x and y as being continuous and restricted to the range — tt to tt. It 
should be remarked that for any fixed value of m, the difference between ^^(m) and 
the ^-function evaluated at r = 2nih/N, will be a nonzero constant that decreases 
with increasing m. 

We can now compare (49) with some results for the stream function on the torus 
derived in earlier times. GreenhillA39 and HicksA40 give the stream function in a 
rectangle. The expression on the torus can be found in the intermediate calculation, 
but perhaps the easiest method of proceeding is the following. 

The stream function on the torus for a vortex at the origin is given as the image 
expression, (cf A40) , 



ij{x,y) 



K 

47r 



E'" 



M 



(a;/27r + Mi)A2 + {y/2Ti + M2)A2 
(MiA2 + M2A2) ' 



where the sums run over all the integers and x and y are restricted to the range — tt 
to TT. Up to an additional constant, 

nx,y) - -l^-Qiy}-. ((^/27r + Mi)A2 + (y/27r + M2)A2)As/2/ ,^0 
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K d 



x/2'K y/2TT 




(s) 



s=0 



K 

2^' 



Z' 



x/2t{ y/27r 




(0). 



Evaluation of the transformation formulaA38 



'k\-s/2V{s/2)Z 
at s = yields 



hi /i2 




(s) = Ti\s/2-lT{l-s/2)Z 





-hi -/i2 



(2-.) 



Z' 



hi /i2 




(0) 



27r 





-hi -/i2 



(2). 



since 



hi /i2 




if /ii and /i2 are not integers. Hence 

ip{x,y) = 



K 



(27r)A2 



(0) = 0, 





a;/27r y/27r 



(2). 



(50) 



Looking at (49), it can be seen that the discrete stream function Wi{m) tends to 
the torus stream function, '^{x, y), as A^ ^ oo for strength k = — 4/tA2. 

HicksA40 gives an expression for t/j^x, y) in terms of simpler ^-functions. This 
is, up to the usual additive constant. 



K 



^{x,y) = -r- 

4n 



— (xA2 - yX2) + In (ei{z/27i, i)ei{zX * /27r, i)) 



(51) 



where z = x ^- iy and 6i{u) is denoted by H{2Ku) in Hicks and Greenhill. Using 
the Jacobi transformation formula for ^-functions, it is easily checked that "0(^7 u) 
is symmetrical in x and y. Our definitions of ^-functions are those of Oberhettinger 
and MagnusA41 where a brief discussion of the motion of vortex systems can also 
be found in chapter 4. 

Incidentally, KroneckerA42 reduced the Epstein ^-function (49) to a form in 
6'-functions, 





hi /i2 



(2) = 27rA2/iiA2-7rln 



0l{ui,Ui)0i{u2,-UJ2) 



(52) 



where rj is the Dedekind function, r]{uj) = q'Al/12 Yli Xoo{l—qX2n). The variables hi 
and /i2 are to be taken in the range from to 1. For the square torus, ui = UJ2X* = i, 
q = eX—Ti and ui = U2X* — {x + iy)/27r. Comparing (52) with (50) and (51), we see 
that Hicks and Greenhill had earlier obtained an equivalent reduction. Kronecker's 
formula has been rediscovered a number of times since. These formulae can be 
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used for numerical evaluation of t/j, (50), although there exists a form in terms of 
incomplete F-functions that is generally more efficient, except for small x + iy. 

For a single vortex, a glance at equation (21) reveals that wXn vanishes for all 
n and the vortex is thus stationary. The energy is given hj E = w{0) = ^ 1/An{k) 
and the enstrophy is ^2 = NX2 — 1. Both these quantities diverge as A^ ^ oo. 

The stream matrix coefficients for a set of vortices is 

i 

where the Hi are the vortex strengths. For two, equal vortices, setting m2 = —mi 
and K2 = Ki = — 47rA2, we get real coefficients wXn = 2 cos (2/c7rn.mi/A^)/An(A;), 
wXO = so that we can place this configuration on the tetrahedron. The enstrophy 
equals 2(A^A2 — 2) while the next invariant is 4(16 — NX2). 

The systems with A^ from 3 to 31, were evolved in time. We present the results 
in Figure 1 for SU(9) as being typical. 

The initial position of the vortices was at mi = (1, 0). In our view no particular 
pattern emerges. 

In any finite scheme, the localisability of individual vortices is lost, the vorticity 
can become redistributed and it is difficult to model the motion of ideal, point 
vortices in this way unless, possibly, A^ is extremely large. However the evolution 
for SU(31), shown in Fig. 2, offers no evidence for such a trend. 

The results for SU(5) were somewhat atypical and are displayed in Fig. 3. The 
vortices appear to be rotating around each other. 

The mode-truncation calculation starting from the same stream function yields 
results of a generally similar nature. 

A number of other initial configurations were also propagated with, again, 
roughly comparable outcomes. 

5. STATISTICAL BEHAVIOUR. 

The numerical results so far presented are for relatively short time evolution. It 
has been suggested by KraichnanA43 that two-dimensional turbulence can be sta- 
tistically modelled on the assumption that the system is ergodic and can, after a 
sufficiently long time be described by a microcanonical or even a canonical dis- 
tribution. In the latter case, two "temperatures" can be introduced as Lagrange 
multipliers for the conserved quadratic quantities, the energy E and the enstrophy 
n. 

A number of computer experiments have been performed in both the viscid and 
inviscid cases on the truncated versions of the Euler and Navier-Stokes equations 
to test this idea. The results are suggestive but not conclusive. 

The models discussed in the present work allow a similar numerical analysis. 
These systems, having more than just the two conserved quantities, E and f2, of 
the truncated theory, might provide a more realistic arena in which to explore the 
statistical hypothesis. 
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With this in mind the systems were evolved for long times, at a reduced numer- 
ical accuracy for speed purposes. The quantities, E and the Si were found typically 
to be conserved to one part in 10A5 to 10A6 over the extended time period. Vor- 
ticity moments and correlation functions were evaluated by simple time averaging 
since actual ensemble averaging was impractical. The results were compared with 
the canonical distribution values, where appropriate. 

The basic theory can be found in Kells and OrszagA44, for example, and so we 
need not give many details. 

The canonical vorticity distribution is given by 

P(C) = ^exp ( - 5^(aAnA-l +/3)|CAn|A2) 



with the partition function, Z = J P{QdC. The relation with the parameters in Ref 
44 is a = CA-1 and (3 = DX-1. 
Ensemble averages are 



(FiO) = J F{OP{C)dC. 



In terms of the two temperatures a and (3 the moments of the vorticity are 

(2r) - IV 
((CAn)Ap) = ^-^^^^ (aA„A-l + /3) X-p (53) 

The time-averaged moments are 



Mp{n,T) = ^ J XT\CXn\Xpdt, 



which, for large T, are to be compared with the ensemble averages (53). 
The mode correlation function C{s) is defined by 

C(n, T; s) = -—— — XT- sCXn(t)CXn(t + s)dt. 

^ ' ' ^ TM2{n,T -s^ ' s V ;s V ; 





The prefactor is a normalisation. If C{s) does not tend to zero with increasing s the 
system is probably not ergodic and cannot be described by a statistical ensemble. 
Three starting distributions were chosen. One was the double vortex discussed 
in section 4, another was a vortex-antivortex pair and a third was a more or less 
random arrangement. In the latter case the coeffients in the truncated model were 
adjusted to give the same energy and enstrophy as the corresponding SU(A^) model 
so that the calculated canonical temperatures a and (3 should be the same. 
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Figs. 4, and 5 display some results using the double vortex starting configuration 
for SU(9). The figure captions are descriptive. Fig. 6 shows the long time evolution 
of the corresponding stream function. For the double vortex the values of the two 
temperatures were calculated to be a = 2.9255188 x lOA-2 and (3 = 0.24223645. 

The presentation is limited to these data sets for reasons of space and also 
because one would really like to explore much larger A^ values. More extensive data 
and discussion is contained in the Manchester PhD thesis of A.W. 

As a measure of the accuracy of the evolution algorithms, we tested the con- 
served quantities. The following are the values of the SU(9) energy and the first 
three Casimirs, 5*2, S3 and 54, at t = (the first number in the brackets) and at t = 
1000 (the second number): ^(14.7549152, 14.7549145), >52 (158.000037, 158.000058), 
>53 (-308.0002, -308.0035), ^4 (38022.0096, 38022.0630). The integer parts of the Si 
are the exact values. The evaluation of the highest Casimirs from the stream func- 
tion coefficients is very time consuming due to the multiple summations. 

The statistical results are inconclusive. We find no indication of non-ergodic 
behaviour but the evidence for a statistical description is still only suggestive. The 
results are perhaps better than one would expect for systems with a small number 
of modes when treated by a canonical distribution which, moreover, ignores the 
other conserved quantities. The micro canonical distribution would be more relevant 
but the evaluation of the statistical averages is then itself an involved numerical 
procedure which we have chosen to avoid. 

The trunction results differ in no essential aspect from the corresponding SU(A^) 
cases except that the vorticity second moment of the highest mode does not relax 
to the canonical distribution value for large times. 

There is also nothing in particular that distinguishes the other starting con- 
figurations although our numerical experiments are not yet very extensive in this 
respect. 

The statistical mechanics of systems of this type, with many conserved quan- 
tities, remains to be elaborated. ZeitlinAlO makes some relevant remarks on the 
structure of the phase space. 
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6. COMMENTS AND CONCLUSION. 

In this paper we have been mainly concerned to set up the SU(A^) models and 
to discuss the nature of the infinite A^ limit, which is not obvious. We have also 
presented a sample selection of numerical results which should be regarded as ex- 
ploratory rather than definitive. 

Our numerical approach was unsophisticated. For the truncated system (which 
is sometimes referred to as the Galerkin approximation) some impressive speed 
advantages could have been achieved by using pseudo-spectralA45 or collocation 
methods combined with fast Fourier transforms. This is, apparently, not possible 
in the SU(A^) models because of the nonlocal coupling terms in the coordinate- 
space form (25). The whole point of the pseudo-spectral method is to calculate the 
mode-coupling terms in the representation in which they are local. For the SU(A^) 
cases we do not seem to have this option. Without some technique corresponding 
to the pseudo-spectral one, the SU(A^) models could never be viewed as numerical 
alternatives to the standard truncation or finite element methods. 

It might be expected that, as A^ becomes bigger, the results for the group 
model and those for the truncated system should approach one another. There 
was no evidence for this in the short-time evolutions that we have performed. Also 
there was no indication that the quantities conserved in the continuum theory were 
progressively better conserved in the truncated versions as A^ increased. Perhaps 
the values of A^ are still too small or it might be that the N ^ oo limit has not 
been closely enough considered and that the expectation is unfounded. 

More disturbing is the oscillatory behaviour of an "entropy" J^n^^sCn/^n as 
a function of time. These evaluations are at a preliminary stage and have not been 
displayed. They may indicate that the systems are not ergodic or that the number 
of modes is small. 

A corresponding analysis can be performed on the two-sphere. Although the 
coupling coefficients are more complicated, the eigenvalues are simpler, being the 
same as in the N —^ oo limit. A discussion will be presented elsewhere. 

Calculations have also been done on a triangular lattice, corresponding to a 
regularly slanted torus. For real mode coefficients, the results are relevant for motion 
on the surface of a regular tetrahedron whereas the square torus discussed in this 
paper gives a degenerately fiat tetrahedron. It would not model the Earth too well. 
This tetrahedron is, in fact one of the fiat Riemann surfaces discussed elsewhereA37 
and it would be possible to extend the present calculation to these. 

Whether or not this whole class of models proves to be of use in approximating 
continuous theories, they at least provide an interesting set of dynamical systems. 

More realistically the effect of viscosity could be investigated by analysing the 
Navier-Stokes equation. 

Another interesting question concerns the Lagrangian stability of the motion 
i.e. of the trajectories in the Lie group. It is known that the Eulerian (or Lie algebra) 
motion can be stable yet the Lagrangian one unstable in the continuous case, being 
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related to the sectional curvature of the group of area-preserving diffeomorphisms. 
It is amusing to consider the discrete analogue of this. 

The computations were carried out on a Hewlett-Packard workstation. Trans- 
ference to a more powerful machine is planned and it is hoped to reach large values 
ofN. 

The recent preprint by RankinA46 contains some information on the SU(oo) 
limit, mostly in a particle physics context. 

After this work was completed we were apprised of the paper by Miller et al 
A47 in which this finite class of models is also discussed. No numerical calculations 
are given and there are a number of formal differences with our setup particularly 
the choice of the discrete Laplacian. 
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FIGURE CAPTIONS. 

Fig.l Double vortex stream function for SU(9) plotted on the torus at evolution 
times t =0, 1, 2, 3, 4 and 5. 

Fig. 2 Double vortex stream function for SU(31) plotted on the torus at evolution 
times t =0, 0.5, 1.0 and 1.5. 

Fig. 3 Double vortex stream function for SU(5) plotted on the torus at evolution 
times t —0, 1, 2, 3, 4 and 5. 
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Fig. 4 SU(9) second moment of vorticity M2 =< CnA2 >t, as a function of energy- 
eigenvalue, An, calculated (by time averaging) at several evolution times, t. The 
diamonds correspond to the initial values, the crosses to t = 750 and the stars 
to t = 1500. The circles indicate the values calculated on the basis of a two- 
temperature canonical distribution. The initial configuration was the two vortex 
one. The relaxation to the canonical values is better than might have been expected 
in view of the relatively small number of modes. 

Fig. 5 Second moment of vorticity and correlation function for three modes of the 
SU(9) model as a function of time {kX2 = nA2). The horizontal levels indicate the 
value expected from a canonical distribution. Note the start of the vertical axis. 
Although the ordering of the levels is not reproduced, the general agreement is more 
than acceptable. The behaviour of the correlation function is consistent with an 
ergodic development. 

Fig. 6 SU(9) long time evolution of the stream function starting from the double 
vortex configuration. The evolution times are t = 200, 400, 600 and 1000. See Fig.l 
for the short time behaviour. 
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